Wavepacket Dynamics in Nonlinear Schrodinger Equations 
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Coherent states play an important role in quantum mechanics because of their unique proper- 
ties under time evolution. Here we explore this concept for one-dimensional repulsive nonlinear 
Schrodinger equations, which describe weakly interacting Bose-Einstein condensates or light prop- 
agation in a nonlinear medium. It is shown that the dynamics of phase-space translations of the 
ground state of a harmonic potential is quite simple: the center follows a classical trajectory whereas 
its shape does not vary in time. The parabolic potential is the only one that satisfies this property. 
We study the time evolution of these nonlinear coherent states under perturbations of their shape, 
or of the confining potential. A rich variety of effects emerges. In particular, in the presence of 
anharmonicities, we observe that the packet splits into two distinct components. A fraction of the 
condensate is transferred towards uncoherent high-energy modes, while the amplitude of oscillation 
of the remaining coherent component is damped towards the bottom of the well. 

PACS numbers: 05.60.Gg;67.85.De;42.50.Md 



I. INTRODUCTION 

Coherent states were introduced in quantum mechan- 
ics by Schrodinger in 1926 to describe minimum uncer- 
tainty wave packets that satisfy the correspondence prin- 
ciple. The standard coherent states are defined as trans- 
lations of the Gaussian ground state of the harmonic 
oscillator potential. The peculiarity of those states is 
that, during the time evolution in such a potential, they 
remain of minimum uncertainty at all times. This re- 
markable quasi-classical evolution is highly non trivial in 
quantum mechanics, the general rule being the spreading 
of the wave packet and the delocalization of the proba- 
bility density. The harmonic oscillator coherent states 
arise in systems whose dynamical symmetry group is the 
Heisenberg-Weyl group. They can be generalized to sys- 
tems with different symmetry groups, like the SU(2) spin 
coherent states, and appear in a wide range of physical 
situations [U [2] ■ 

If an initial Gaussian wave packet is subjected to the 
action of an anharmonic potential, it will generally spread 
out. In some cases, after the initial spreading, the quan- 
tum state may, periodically, come back almost com- 
pletely to its initial state. This revival of the wave packet 
occurs in systems where the spectrum may be expanded 
locally in terms of a quantum number, a characteristic 
situation of one-dimensional integrable Hamiltonian sys- 
tems [21|3]. In contrast, if the corresponding classical dy- 
namics is chaotic, the wave packet will spread and relax 
towards the phase space chaotic region, with time depen- 
dent fluctuations of the density that reflect interference 
effects. The structure of the underlying classical Hamil- 
tonian thus has a strong influence on the dynamics of the 
packet, and may produce quite different effects depend- 
ing on the integrable or chaotic nature of the classical 



dynamics [3]. 

Here we are interested in a situation where the clas- 
sical dynamics is simple, we consider integrable one- 
dimensional Hamiltonian systems. However, the diffi- 
culty is related to the more general character of the 
quantum dynamics considered, since we include nonlin- 
ear terms in the Schrodinger equation. The resulting non- 
linear Schrodinger equation (the Gross-Pitaevskii equa- 
tion, GPE) has a wide range of physical applications. It 
emerges, in particular, in two important cases: in the de- 
scription of a Bose-Einstein condensate (BEG) of weakly 
interacting particles [5| , and in the description of electro- 
magnetic waves (light) propagating through a nonlinear 
medium [B]. 

The first point we are interested in is to determine if, in 
the nonlinear case, there still exist coherent states, in the 
sense of a set of initial states that are able to propagate 
in time without spreading or changing their shape. This 
question is particularly relevant in the context of BEGs, 
since the mere existence of a coherent motion means, 
physically, that the condensate is preserved in time and 
that the atoms do not diffuse to different modes during 
the motion. We consider here the particular case of a pos- 
itive nonlinear coefficient, which corresponds to a BEG of 
repulsive interactions, or to a defocusing medium in non- 
linear optics. The most elementary expectation would 
be that the additional repulsive nonlinear term in the 
Schrodinger equation enhances the spreading of an initial 
wave packet. This is of course true for the free propaga- 
tion. However, as in the case of the linear Schrodinger 
equation, we find that a particular role is played by the 
harmonic confining potential. For that potential it is 
shown that the phase-space translations of the nonlinear 
ground state behave as coherent states, e.g. during the 
time evolution the center of the packet follows a classical 



phase-space trajectory, without any change of its shape. 
These translations therefore constitute a set of nonhn- 
ear coherent states which will be properly defined in sec- 
tioning This behavior is specific of the harmonic poten- 
tial. Furthermore, we study the stability of the nonlinear 
coherent states under deformations of their shape. For 
small deformations, the packet remains coherent and its 
center follows the corresponding classical trajectory, with 
superimposed small shape oscillations of frequency given 
by the multipole modes of the ground state. Remarkably, 
this result holds also for large initial perturbations. For 
instance, the motion of a very compressed initial Gaus- 
sian state can be decomposed into a standard dipolar 
motion of its center and a superimposed large amplitude 
shape expansion and compression cycle. 

The next relevant question concerns the evolution of 
a nonlinear coherent state subjected to an arbitrary ID 
confining potential. In contrast to the linear case, when 
both anharmonicities and nonlinearities are present the 
spreading and revival of the packet are not observed, and 
a new phenomenology emerges. We find, as in previous 
studies [7], that for small anharmonicities and small am- 
plitudes of oscillation the packet keeps, to a good ap- 
proximation, its coherence. Its center follows a classical 
trajectory with superimposed small shape fluctuations. 
However, as the anharmonicity or the amplitude increase, 
a new process appears. The packet splits into two com- 
ponents, where part of the packet is damped towards the 
bottom of the potential, while the rest leaves the packet 
to form an uncoherent higher energy phase-space cloud. 



II. WAVEPACKET DYNAMICS IN A 
HARMONIC POTENTIAL 

A. Coherent states of the Gross-Pitaevskii 
equation 

We consider the one dimensional time-dependent 
Gross-Pitaevskii equation: 






h^ a2^(x,t) 



V{x)^'{x,t) 



2m dx^ 

(giV|vE'(a;,t)|2-^)*(x,i) 



(1) 



which describes, in the mean-field approximation, the 
dynamics of a Bose-Einstein condensate of N identical 
bosons, in the presence of repulsive interactions (5 > 0), 
in an external potential V{x) [5 . Here, '^{x,t) is the 
normalized wavefunction of the condensate, m the mass 
of each particle, g the interaction constant, and jj. the 
chemical potential. Aside from cold atom physics, it 
has been shown that Eq. M provides an accurate de- 
scription of many interesting physical problems, among 
which we can mention hydrodynamics [5], or nonlinear 
optics 6]. In the latter case, the ID GPE can be derived 
from the propagation of light in a two-dimensional non- 
linear medium, under both the monochromatic and the 
paraxial approximations. 



We assume that V{\x\ — >■ 00) — >■ 00, and look for solu- 
tions of the GPE which evolve in time without changing 
their shape. We thus seek for solutions in the form 



^(x, t) = (j){x—xt){t), t) exp 



ipo{t) 



xa{t) 



, (2) 



where xo{t) and poit) are real functions of time. This so- 
lution represents a time dependent evolution in which the 
wavefunction is translated along the phase-space trajec- 
tory {xo{t),P(){t)). The substitution of Eq.([2| into Eq.Q 



gives 



ih 



dt 



X — XQ{t),t 



2m dx^ 



+ V{x)<P{x~xo{t),t) 

X~XQ{t),t 

+ {9N\(l){x - xo{t),t)f - fi) (p{x - xo{t), t) 

Po{t)\ d(j)_ 
dx 



+ ih xo{t) 



2 \ m 






+ Po(0 2; 



Xo{t) 



{x-XQ{t),t) 



(3) 



where Xoit) = dxo/dt and po(t) = dpo/dt. Equation 
([3]) takes a simpler form if the phase space trajectory 
(xo{t),Poit)) coincides with a trajectory of the corre- 
sponding classical non-interacting problem 



Xoit) = 

Po{t) = 



Pojt) 
m 

dx 



Xo{t) 



Making the change of notation x — xo{t) — )■ x, Eq. (Is]) 
simplifies to 






V{x + xoit))- 



dV 
dx 



^o(t) 



Xo{t)^ 



(4) 



in which (j) and its derivatives are now evaluated in (x, t). 
Equation Q shows that, in the reference frame of the 
classical trajectory, the particle feels a time- dependent 
potential. In the new reference frame, the coherent state 
should be a stationary state of Eq. Q. The stationarity 
condition imposes a time-independent potential. This 
leads, for any x, to the condition 



d 
dxo 



V{x + xo) 



dV 
dx 



(-?) 



In particular, for 2; = it takes the form 



d^V dV 



Xo 



dxl 



dxo 



= 



(5) 



(6) 



This equation is satisfied if and only if V{x) is a quadratic 
function of x. Hence, the only function that produces, in 
the new reference frame, a time-independent potential is 
the harmonic one. Finally, for a harmonic potential and 
in the reference frame that follows the classical phase- 
space trajectory, the quantum equation of motion takes 
the form 



ih 



dt 



2m dx' 



+ V{x)cb+{gN\cl>\^~ti)<j>. (7) 



Therefore, the coherent states of the GPE are defined by 
its stationary states, that satisfy the equation 

~^^ + ^^'^^'f'o + 9N\M'=^^^o■ (8) 
It follows that, for a harmonic potential. 



*o(a;, t) = (f>o{x - xo{t)) exp 



fipo{t) 



xo{t) 



\ h \ 2 

■ ■ (9) 

is a time-dependent exact solution of Eq. (fl). Here, 
{xo{t),po{t)) is a phase-space trajectory of the corre- 
sponding non interacting classical system. In other 
words, the time evolution of the wave packet defined by 
Eq. ^ reduces simply to the time evolution of its center, 
that follows a classical trajectory. Among the different 
possible stationary states ^o of Eq. ([s]), it is customary 
to define as the standard coherent state the ground state, 
which minimizes the energy as well as its spatial exten- 
sion ^. From now on we refer to the set "^olx^t), with 
(/>o defined as the ground-state of Eq. (Is]) and xo(0) and 
Po{0) arbitrary, as the set of nonlinear coherent states. 

It is easy to see that the previous results are not 
only valid for a quadratic nonlinearity of the GPE, 
but that they hold in fact for an arbitrary exponent 
~ gN\'i{x,t)\"-. This remark extends our results to a 
large family of nonlinear Schrodinger equations. 

In order to illustrate the previous results, we have nu- 
merically computed the time evolution of Eq. (fTl), and 
plotted the phase-space Husimi distribution of the wave- 
function at different times. This distribution is defined 
as 

Hix,p,t)^\{xp\^m^ , 

where \x p) is a standard linear harmonic oscillator co- 
herent state centered around the phase-space point {x,p), 
whose X representation reads: 



{x\xo po 



>-{ 



mio 



nh) 



1/4 



exp 



{x- xof 



X exp 



«Po 



(^-t) 



(10) 



The typical width of a standard coherent state in the 
X and p directions is Xho = {2h/mujy'^ and pno = 
{2hmuj) ' , respectively. To obtain Fig. ([ij), we nu- 
merically calculate the ground state (jioi^) of the Gross- 
Pitaevskii equation in a harmonic trap, V{x) = ^muj'^x'^ , 
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FIG. 1: (color online) Time evolution of a shifted ground 
state of the GPE with a harmonic confining potential, with 
parameters 7 = 115, d = 5. Husimi representations of the 
wavefunction are given at times ujt/{2-n) — (a), 0.25 (b), 0.5 
(c), 0.75 (d). The (red) full curve is the classical trajectory 
of the corresponding linear problem of energy given by the 
center of the initial packet. 



and then compute the time evolution of a translated 
ground state, ^(x, i = 0) = (f)o{x + d). In order to charac- 
terize the intensity of the nonlinearity, it is convenient to 
define a dimensionless parameter. In terms of the char- 
acteristic width Xho and energy huj of the ground state 
of the noninteracting harmonic oscillator, we define the 
parameter 7 = 2gN/(x-aof^'^), 



7^ 



12m gN 



(11) 



As predicted above, in the nonlinear case the wave 
packet dynamics reduces to a simple phase-space transla- 
tion of its center, that follows the corresponding classical 
trajectory (full (red) curve in the figure). During this 
process, its shape does not vary in time, and there is no 
rotation either. In particular, the shape of its projection 
onto the x axis does not change in time. 

This behavior qualitatively differs from the dynam- 
ics of the linear Schrodinger equation (non-interacting 
case), where the motion of an arbitrary initial wave- 
function in a harmonic trap consists in a phase-space 
rigid rotation with respect to the origin |10j : defining 
z — x/xno + ip/pnoj it is known that the linear evolu- 
tion of an arbitrary initial Husimi distribution Hoiz) in 
a harmonic oscillator reads: 



niz,t) = nix,p,t)=naize''^'). 



(12) 



This implies a rigid phase-space rotation of any initial 
state. 
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FIG. 2: (color online) Time evolution of the linear Schrodinger 
equation with a harmonic potential with parameters 7 = 0, 
d = 5. The initial state is the same one as in Fig. ([T| (a 
shifted groundstate of the GPE) .(a): Husimi representations 
of the wavefunction are given at times a;i/(27r) = (a), 0.25 
(b), 0.5 (c), 0.75 (d). 



To stress the difference between the linear and the non- 
linear dynamics, we plot in Fig. ([2| the linear evolution of 
the same initial state as in Fig. TO. We observe that, in 
contrast to the nonlinear evolution, the initial packet now 
rotates as it follows the classical trajectory, and therefore 
changes its shape as a function of time in the position 
representation. The coherent state of the linear case cor- 
responds, necessarily, to a perfectly spherical Gaussian 
initial packet, a shape which is invariant under rotations 
in any representation. 

Remark that the classical trajectory followed by the 
center of the packet has no dependence on the interaction 
parameter g. It is a classical trajectory of the noninter- 
acting problem, fixed by the initial position of the packet. 
In particular, the frequency of the oscillation is indepen- 
dent of the interaction, a result demonstrated by Kohn 
[llj for the cyclotron frequency of interacting particles, 
that was later on generalized to interacting particles in a 
parabolic confining potential |12j . 

The experimental realization of nonlinear coherent 
states, as well as the control of their initial phase-space 
location, is a natural procedure in the context of cold 
atom physics. This is because cold atoms are usually 
trapped in parabolic magnetic potentials, and the cor- 
responding BEG is thus a coherent state centered at the 
bottom of the potential. Phase-space translations of that 
state are easily implemented by a sudden shift of the 
trap with respect to the condensate. The study of dipo- 
lar oscillations were among the first experimental tests 
of excited collective states [T3] • More recently, dipole ex- 
citations were used to test transport properties of BEGs 



across an impurity [TIHTB] or through disordered poten- 
tials [m Uni [17] . Dipole oscillations were also proposed 
as a test of the existence of a superfluid phase for light 
moving in a nonlinear medium |18| . 

The quantum dynamics in the presence of nonlineari- 
ties is thus particularly simple if the initial state is a co- 
herent state. What happens to an arbitrary initial state? 
We will explore in detail this question in the following sec- 
tions, which will be particularly relevant in the context 
of nonlinear optics since, in contrast to BEGs, in optics 
gaussians are the natural transverse intensity profiles. 



B. Stability of the oscillations 

In this section, we study the stability under defor- 
mations of the initial wave packet '^o{x,t) (Eq. ([9])) in 
the presence of a harmonic confining potential V[x) — 
^rauP'x'^ . For that purpose, we look for solutions of the 
GPE having the form (l2| and were (l){x,t) = 4>o{x) + 
(5(/)(a;, t). Actually, the problem of the stability of the time 
dependent solution ^0(2;, i^ equivalent to the problem 
of stability of the stationary ground state of Eq. (fTl) . The 
first order expansion in 5(j) of Eq. (It]) leads to 



ih 



dS(j) 

Ik 



2m dx"^ 



V{x)5cj) - fiS(f) 



-gN{2\cj)ofd^ + cl>lS^ 



(13) 



which, with its complex conjugate equation, form the so 
called Bogoliubov-de Gennes (BdG) system. Since 0o is 
real, the BdG system reduces to 



ih 



d_ 
dt 



= M 



(14) 



where M 



and A 



A gN4,l 
-gN4>l -A 

The stability of the solution 



hf_jr_ 

2m dx^ 



-V(x) 



t>Q IS given 



by the sign of the eigenvalues hujn of M, which are the 
energies of the elementary excitations [ UmVn 1 , given by 



huJn 



^ M 



(15) 



Our calculations are the ID equivalent of the 2D work 
of Ref. [H], and we will not give the technical details 
here. For instance, in the strongly interacting limit (the 
so called Thomas- Fermi limit) , the spectrum is given, for 

n e N* , by 



LU 



n{n + 1) 



(16) 



This result shows that the frequencies become, in that 
limit, independent of the nonlinearity and that the n = \ 
dipolar excitation is unchanged, uji = u. All eigenvalues 
are real, a fact that ensures the dynamical stability of 
the coherent state under small deformations. 



In the following we use a different method to test the 
stability of the motion of coherent states under shape 
deformations. We use the Virial theorem for the GPE 
[201121] and, applying a variational principle, recover the 
former results as well as some extensions of their regime 
of validity. The Virial Theorem states that, for a solution 
^(a;,i) of Eq.Q, the average spatial extension < x^ > 
of 'i(x,t) verifies 



dt < x' 



>= 



4E 



K 



"^EpfL — 2 



dV 
dx 



where 



Ek = / — |9,*(x,t)pdx, 
J 2m 

Enl = ^ J \^{x,t)\Ux , 
Ep = f V{x)\^{x,t)\^dx , 



(17) 

(18) 
(19) 
(20) 



and dx 



dt 



< A{x) >= J A{x)\^ix,t)\^dx 



dx] '-'t — dt' ^ 

for any function A{x). This theorem has been used in 
particular to study the collapse dynamics of a BEG. It 
is important to mention that Eq. (17) follows from the 
fact that ^ extremizes the Gross-Pitaevskii functional 
Em = Ek + Enl + Ep. The quantity E = Ek + 
Enl + Ep does not depend on time. In the particular 
case V{x) 
to 



2^2 



jmoj X 



the relation (x^) = 2Ep leads 



d' < x' 



>= - [4Ek 
m 



2E 



NL 



AEf 



(21) 



For instance, for the non-interacting case, 5 = 0, Ej^l ~ 
0, and thus E = Ek + Ep is a constant determined by 



the initial condition. Then, Eq.(21) simplifies to 



d'i < x^ >-- 



-Auj^ 



<x"- > 



E 



(22) 



This means that, for any initial wavefunction, the spatial 
extension of ^(a;,i) is an oscillatory function of time, 
with frequency 2a;, a fact clearly seen in Fig. ([2]). Indeed, 
since, as we mentioned previously, the dynamics in the 
non-interacting (linear) case of a harmonic oscillator is 
simply a rigid rotation in phase space, it is clear that 
every half-period of the oscillator the spatial extension 
comes back to its initial value. 

We now take into account the presence of interactions, 
and more particularly, we assume to be in the Thomas- 
Fermi limit 7 — > oo. The reason for this assumption 
is that in this case an explicit form of the ground state 
(j>Q^{x) is known 



0r(^) 



I fi — ^mui'^x'^ 



gN 



(23) 



(for x^ ^ 2/i/mu'^, whereas (j)^^{x) — for x'^ > 
2ii/mw'^). In order to solve Eq. (21 1, we assume that 



the wavefunction is able, during its time evolution, to 
follow the classical trajectory as well as to vary its spa- 
tial extension, denoted L. For |a; — a;o(i)| < L{t) we write 
it in the form 



^!Lix,t)^C{L{t)) 



X exp 




(24) 



and ^L(a;,t) = if |a; - xo{t)\ > L{t). In the latter 
expression, C{L) = y/S/iL ensures the normalization of 
'$L{x,t) at any time. Let us substitute Eq.(24) into the 



virial theorem (21), in which all terms depend only on 
L(t), Xo(t) and Po(t), respectively noted L, Xq and po for 
a matter of readability, and their derivatives xq = '"S^ , 
9psM and L = ^^ 



Po 



dt 



dt 



2 

, pI 



L 



^2^ 



2{xqXo + Xq^ 



L' 



-2u\xi + —) + 



3gN 



(25) 



Using the classical equations of motion, all the terms con- 
taining information concerning the classical trajectory 
vanish, and we finally obtain 



LL 



L^ = -Lj^L^ 



3g7V 
2mL 



(26) 



The equilibrium solution of the latter differential equa- 

1/3 



tion is Lp 



\ 2mu:^ J 



which coincides with the usual 



^eq — y2muj^ 

spatial extension of the Thomas-Fermi solution. Let us 
now consider small deviations with respect to its ex- 



tension, and write L(t) in Eq.(26) in the form L{t) 
Leg + 5L(t). Performing a first order expansion in u{t) 
SL{t)/L^g < 1, we get 



3uj u 







(27) 



which describes a periodic oscillatory motion of the width 
of the wave packet of frequency y/iuj. This frequency cor- 
responds to the n = 2 quadrupole mode of the excitation 



spectrum of Eq.(16). To summarize, in the two limiting 
situations 7 = and 7 -^ cx) the quadrupole deforma- 
tions of the time dependent coherent state are stable and 
the corresponding frequencies are 2a; and •s/S'^, respec- 
tively. 

In order to study the intermediate regime, for which 
we have no analytical expression of the ground state, 
we choose to use a normalized gaussian ansatz $''(x,t) 
(which tends to the correct form in the absence of non- 
linearities), with a time-dependent width 77 (i) 



$"(0;,*) 



1 



(2^7y2(f))i/4 

'ipo{t) 
X exp ' 



exp 



{x-xo{t)f 

Xo{t) 



(28) 



The same procedure as before leads to the following dif- 
ferential equation for 77 (i) 



2{T]ii + 77 



2^2 



2m?T] 



2„2 



2a;^?7 



gN 



2 ^TT 771 77 



(29) 

Replacing in Eq.(28l the stationary width 779 = a^Ho/2 
of the linear g — Q limit of Eq.(29) gives the function 



<^'^°{x,t) which coincides with the well known definition 
of the usual coherent state of the harmonic oscillator, de- 
fined by the complex parameter z — xq/x-hq -I- ipo/pao- 
For a non zero interaction constant, u{t) = ri{t)/r]Q veri- 
fies 



uii + u = u 



1 



^1 



(30) 



Let us denote ^£^(7) the strictly positive equilibrium so- 



lution of Eq.(30|. Ueqij) is an increasing function of 7, 



equal to 1 for 7 = 0, and that tends to infinity in the 
limit 7 — )■ 00. Similarly as above, we perform a first 
order expansion writing u{t) = u^q{j) + Su(t), and as- 
suming Su{t) ^ Ueqij), to obtain again a second order 
differential equation 



Si 



n^du = 



(31) 



where Q is, in this approximation, the quadrupole fre- 
quency 



n^ = i 



1 



Uegh)^ 



(32) 



Note that equations ( 30 ) and ( 32 ) have been already ob- 



tained by a variational principle in Ref. [7] including also 
the fourth order moment as time-dependent parameter. 
In the linear limit 7 = 0, Ueq{0) = 1, and we recover 
ri = 2w, as it should. In the other limit of strong non- 
linearity, Ueqij -^ 00) -^ 00, and we recover il = a;v3, 
which is the correct result, as was shown previously. In 
Fig. ([3]) we plot the comparison of Eq. ( 32 ) for arbitrary 7 
to a numerical calculation of the quadrupole frequency. 
Despite the fact that the gaussian ansatz is only cor- 
rect in the linear limit, we see that it provides a quite 
good approximation of the quadrupole frequency for ar- 
bitrary 7. (The numerical simulation performed is the 
following: For any 7, we numerically compute the ground 
state of the GPE, we spatially shift it from the bottom 
of the potential, and apply a (norm-preserving) deforma- 
tion (getting SL{t = 0) = O.lLeq) in order to excite the 
quadrupole mode. Then, the real-time evolution of the 
GPE is computed, and the frequency which maximizes 
the Fourier transform on SL{t) is finally found.) 

The previous results show the stability of the coher- 
ent states (and therefore of a condensate) under small 
shape perturbations when moving in a harmonic poten- 
tial, and provide the typical frequencies involved. We 
have also numerically explored the evolution of packets 
whose initial shape strongly deviates from the coherent 
state. For instance, in Fig. Q we show the nonlinear 
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FIG. 3: (color online) Normalized square quadrupole fre- 
quency for different values of the nonlinear parameter 7. 
(Black) dashed line represents the numerically computed fre- 
quency, full (blue) line the analytical result obtained using a 
gaussian ansatz variational principle. 



evolution of a Gaussian coherent state of the linear prob- 
lem (defined as the translated Gaussian ground state of 
that problem). What is observed is the usual dipole os- 
cillation following the corresponding classical trajectory 
with a superimposed large amplitude quadrupole vibra- 
tion. The spatial width of the initial Gaussian state is 
small compared to the corresponding nonlinear state, see 
Fig. (IT]). It follows that, because of the repulsive inter- 
actions, the packet strongly spreads in phase-space, pre- 
dominantly in the p direction (particles accelerate, see 
part (b) of the figure). This acceleration produces a spa- 
tial spreading of the packet, whose barycenter follows 
the corresponding classical trajectory (part (c)). At this 
point the expansion stops, compensated by the harmonic 
confinement, and a compression phase follows, to recover 
its initial shape. The process can start again. We have 
numerically computed the period of the expansion and 
compression cycle, and found a period (normalized to the 
harmonic oscillator period) Tg/T ~ 0.551, which is close, 
but nevertheless different, from the quadrupole frequency 
predicted from Fig. ^ for the corresponding value of 7, 
T4/T ~ 0.575. 



III. ANHARMONIC EXTERNAL POTENTIAL 

We now explore the robustness of the motion of nonlin- 
ear coherent states when the considered potential differs 
from the harmonic oscillator. More generally, we wish 
to explore the nonlinear motion of initial wave packets 
under an arbitrary potential. Experimentally this is a 
relevant problem since anharmonic potentials are either 
use on purpose [H], or they come as corrections to the 
nearlly harmonic usual traps. From a theoretical point 
of view, frequency shifts and coupling of collective modes 
due to anharmonicities was explicitly investigated in the 
past [7]. 
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FIG. 4: (color online) Time evolution of the GPE with a harmonic confining potential with parameter 7 = 115. The initial 
state is the shifted Gaussian ground state of the linear problem. Husimi representations of the wavefunction are given at times 
u}t/{2-K) = (a), 0.02 (b), 0.13 (c), 0.25 (d), 0.50 (e), 0.51 (f), 0.53 (g), 0.75 (h). 



As an example we consider a potential of the form 



V{x) = -muj a; (1 + ax ) 



(33) 



where a controls the strength of the anharmonicity. We 
consider as initial state the nonlinear coherent state of 
the corresponding harmonic oscillator, i.e. we compute 
the ground state of the nonlinear equation with a = 
(the use of the true ground state does not qualitatively 
modify the results). This state is then shifted along the 
X direction in order to locate the center of the packet 
at X = —d, with d positive. The time evolution of such 
state is then computed for the full potential including the 
quartic term. As d increases, the strength of the quartic 
term of the potential compared to the harmonic one in- 
creases. This strength is measured by the dimensionless 
parameter /3 = ad^ . We thus study how the dynamics of 
the initial packet changes as a function of /3. 

Figure (Isl) shows the time evolution for /3 — 0.04. Be- 
fore analyzing the results, it is useful to show the time 
evolution in the linear case. In the absence of nonlinear 
terms in the Schrodinger equation the time evolution is 
made of cycles of spreadings of the wave packet followed 
by a revival, i.e. after the spreading the packet comes 
back, to a good approximation, to its initial state, and 
the process starts again. This is indeed what is observed 
when 7 = for an arbitrary value of /3, see upper part of 
Fig. (§. 

The motion of the corresponding coherent state in the 
presence of nonlinearities is quite different. For small 
values of /3, such as the bottom panel of Fig.(l5]), we ob- 
serve that the nonlinear dynamics is more robust than 



the linear one. For such values of f3 no spreading is ob- 
served. The packet keeps, to a good approximation, its 
initial shape during the time evolution, while the center 
follows the classical trajectory. Small amplitude dipole 
oscillations are observed, as well as a periodic motion of 
the tilting angle of the axis of the packet with respect to 
the z-axis. But roughly the packet (e.g. the condensate) 
preserves its coherence. 

Things change qualitatively as /3 increases, as shown in 
Fig. ^. For larger initial amplitudes of the oscillation, 
at fixed a, a strong deformation of the packet is observed 
during its time evolution. The packet does not preserve 
anymore its coherence. As it evolves, a filamentary struc- 
ture develops from the packet and winds in the clockwise 
direction around it. This filament extends up to very high 
energies (see part (b) of the figure). By energy and mass 
conservation, the remaining packet has a smallersize and 
its center now occupies classical orbits of smaller energy, 
e.g., its amplitude of oscillation decreases. As time goes 
on, the winding filament compresses towards the packet. 
In this process, the different loops of the filament start 
to interfere. Finally, the reduced packet is completely 
damped at the bottom of the well, and coexists with a 
low density component which occupies a large fraction of 
the higher energy phase space, as shown in part (d) of 
Fig. g. 

This is a remarkable process, that completely differs 
from what is known from the time evolution of the lin- 
ear Schrodinger equation. Using the language of Bose- 
Einstein condensates, one can summarize it as follows (a 
similar effect is expected for, e.g., light motion in a non- 
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FIG. 5: (color online) Time evolution of the linear Schrodinger equation (7 = 0, top panel), and of the GPE (7 = 115, bottom 
panel) in the presence of an anharmonic (quartic) confining potential with parameters a = 0.01, and /3 = 0.04. For each panel, 
the initial state is the corresponding (linear or nonlinear) coherent state, computed for a = 0. Husimi representations of the 
wavefunctions are given at times u)t/(2iT) — (a) and (e), 3 (b) and (f), 12 (c) and (g), 40.5 (d) and (h). 
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FIG. 6: (color online) Time evolution of the GPE in the pres- 
ence of an anharmonic (quartic) confining potential, with pa- 
rameters 7 = 115, a — 0.01, and /3 = 0.5 . Husimi represen- 
tation of the wavefunction at times a;f/(27r) — (a), 1.5 (b), 
4 (c), 40 (d). 




FIG. 7: (color online) The fluidity factor, defined by the ratio 
of the average amplitude of oscillation in a stationary regime 
and the initial amplitude d, versus /3 = atf , for 7 = 115, and 
a — 0.01. The initial state is a shifted nonlinear coherent 
state. The inset shows the same plot for different values of a, 
using for both axis the same scale as the main figure. 



linear medium). In the presence of anharmonicities, the 
kinematic energy stored as center of mass motion of the 
condensate is not preserved, as for a harmonic potential. 
Instead, during the dynamical evolution, one observes the 
emergence of two components. The initial packet is not 



totally destroyed. In the course of time, it looses part of 
its mass, and its amplitude of oscillation diminishes, to 
eventually be almost stopped at the bottom of the po- 
tential. The fraction of the condensate that leaves the 
packet occupies, in the course of time, high energy tra- 
jectories in an uncoherent way. One may speak of some 
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sort of evaporative process, where the initial kinematic 
energy of the center of mass is transformed during the 
time evolution into uncoherent motion of high energy 
particles (evaporation), whereas the remaining fraction 
of the condensate cools down towards the bottom of the 
potential (damping effect). 

At a given evolution time, the amplitude of oscillation 
of the remaining packet depends on /3. To illustrate this 
point, we have computed, as a function of /3, the time av- 
erage of 2[{x){t)/dY, where {x){t) = /^ a;|V'(a;, t)pdx. 
The time average is computed for long times, starting 
from a time such that the evolution of {x){t) looks sta- 
tionary in time. This factor, that we call the fluidity 
factor, is equal to one at /3 = (no damping of the wave 
packet) and equal to zero for a packet totally damped, 
almost at rest at the bottom of the potential. The result 
is represented in Fig. [7J A strong decrease is observed 
as /3 increases. For small values of /3, there is no plateau 
where strictly no damping is observed. The fluctuations 
are due to the interactions between the low-density high- 
energy component with the main wave packet compo- 
nent. It may well be that if we further increase in time 
the position of the time average window the fluidity fac- 
tor globally decreases. That would mean that at very 
long times the packet is always fully damped. We cannot 
give for the moment a definite answer to this point. 

We have also explored the dependence of this process 
on the different parameters. The inset of Fig.JTJshows the 
dependence of the fluidity factor on /3 for packets that 
propagate in potentials with different values of a. The 
superimposition of the curves shows that, on average, 
this quantity depends on a and d only through /3 = ad^. 



IV. CONCLUDING REMARKS 

We have shown the existence of non-spreading states 
for the repulsive GPE, the so called nonlinear coherent 
states. They are defined as phase-space translations of 
the ground state of the nonlinear equation in presence 
of a harmonic confining potential. Due to the repulsive 
interaction, they are strongly elongated in the spatial di- 
rection. In the presence of a harmonic potential, the 
nonlinear coherent states do not vary their shape dur- 
ing the time evolution, their center simply follows a cor- 
responding classical trajectory (of the linear problem). 
This means that the center of mass motion is decoupled 
from other modes of the system. In particular, they are 
stable under shape deformations. We have computed the 
corresponding frequencies of oscillation for different non- 
linearities. In the presence of a harmonic potential, the 
nonlinear coherent states thus preserve their coherence 
during the time evolution. 

The physics is quite different when the nonlinear coher- 
ent states evolve in an anharmonic potential. We found 
that the time evolution now leads to a partial destruction 
of the initial packet (or of the condensate in BECs) . Dur- 
ing the time evolution, the system splits into two compo- 



nents. A fraction of the initial density leaves the packet, 
to occupy high energy phase space trajectories (evapora- 
tive process). The remaining fraction of the packet con- 
tinues to oscillate around the bottom of the well but, by 
energy conservation, its amplitude now decreases (damp- 
ing process). The anharmonicity of the potential thus 
induces a coupling between the dipole mode and other 
excitation modes. The initial center of mass kinematic 
energy is now partially transferred to a fraction of the 
particles, that leave the system, while the amplitude of 
the collective dipole motion of the remaining coherent 
component is damped. This process depends on the an- 
harmonicity and on the initial amplitude through the pa- 
rameter /3, with a stronger damping for stronger values 
of/3. 

In the presence of interactions, the revival phenomenon 
that occurs in linear quantum mechanics thus disappears 
and is replaced by a totally different mechanism. In the 
language of cold atom physics, the condensate is partially 
destroyed and damped when it evolves in an anharmonic 
confining potential. 

Coherent transport and superfluidity are often tested 
by adding an external perturbation, like for instance the 
study of the damping of dipolar oscillations in BECs in 
the presence of an obstacle [TH[Tn]. In ID, the dissipa- 
tive mechanism that breaks superfluidity is related to the 
emission of solitons. Also, loss of coherence and damping 
of collective excitations are predicted as temperature in- 
creases [131 ■ Here, we have shown loss of coherence and 
dissipative effects in the absence of obstacles, simply in- 
duced by the presence of anharmonicities in the confining 
potential. 

Many interesting problems remain open, like a study of 
the motion of initial nonlinear packets in higher dimen- 
sional potentials, integrable or chaotic. The nature of the 
evaporative process described here should be further in- 
vestigated using methods that go beyond the mean field 
approximation. Many studies already exist for the prop- 
agation of ID packets in the presence of random poten- 
tials [T31I21H25]. However, the present one-dimensional 
results show already the deep differences that exist be- 
tween the linear and nonlinear cases in the presence of 
simple potentials. Experimental tests of these differences 
are relatively easy, in particular in the cold atom context, 
by shifting a BEC with respect to an anharmonic poten- 
tial. In this article we have also explored the nonlinear 
dynamics of gaussian wave packets in both harmonic and 
anharmonic potentials, a problem that is relevant in op- 
tics experiments. 
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